Observation of Polarization in Bottomonium Production at ^fs = 38.8 GeV 
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We present a measurement of the polarization observed for bottomonium states produced in p-Cu 
collisions at y/s = 38.8 GeV. The angular distribution of the decay dimuons of the T(1S) state show 
no polarization at small xp and pt but significant positive transverse production polarization for 
either pr > 1.8 GeV/c or for xf > 0.35. The T(2S+3S) (unresolved) states show a large transverse 
production polarization at all values of xf and pr measured. These observations are compared with 
an NRQCD calculation that predicts a transverse polarization in bottomonium production arising 
from quark-antiquark fusion and gluon-gluon fusion diagrams. 

PACS number: 13.88.+e, 14.40.Nd 



It has been known for some time that the observed 
production rates of charmonium and bottomonium res- 
onances in hadronic collisions are much larger than 
the predictions of lowest order Perturbative Quantum- 
Chromodynamics (PQCD) Q. A calculational approach 
based upon Non-Relativistic Quantum Chromodynamics 
(NRQCD) has emerged as a reliable framework for cal- 
culating onium production || . 

Data on the direct production of the charmonium 
mesons ^(1S) and ip(2S) at high energies, when com- 
pared with the predictions of NRQCD, indicate that color 
octet contributions dominate the cross section and that 
S state charmonia are produced through gluon fragmen- 
tation into a 3 S[ 8 ^ octet state ||. Recent investigations 
have shown that the contribution of color octet states to 
onium production may also be very important at fixed 
target energies, but quantitatively the picture is far from 
complete In particular, NRQCD predictions dis- 

agree with measurements of the polarization of ^(1S) and 
-0(2S) mesons produced at collider M and fixed target 
energies |J. 

In NRQCD, the predicted spin effects in onium produc- 
tion can provide further tests of and constraints on the 
various color octet contributions. The quark-antiquark 
fusion and gluon-gluon fusion diagrams which are ex- 



pected to dominate onium production at fixed target 
energies yield significant transverse polarization jjj for 
the produced bottomonium mesons T(1S), T(2S), and 
T(3S). The polarization results in a 1 + cvcos 2 (6*) de- 
cay angle distribution for the polar angle of the decay 
dimuons in the Collins-Soper frame Q. Transversely, 
longitudinally, and unpolarized states decay with a = 
+1, —1, and respectively. 

We have studied the production of dimuons in the col- 
lision of 800 GeV/c protons with a copper beam dump, 



P 



Cu 



n + n + X. 



The apparatus was originally constructed for Experiment 
605 |p and was located in the Meson East Laboratory at 
Fermilab. The data reported here were taken as part of 
a subsequent experiment, Experiment 866 [fiof . Details 
of the apparatus used in E866 and a full description of a 
similar study of the polarization of dimuons from char- 
monium states can be found in Reference [[|. 

Here we present polarizations derived from the an- 
gular distribution of 2 million dimuons in the range 
8.1 < 171^+^- < 15.0 GeV. The data, after analysis 
cuts, cover the kinematic range 0.0 < xp < 0.6 {xp is 
the fractional longitudinal momentum of the dimuon in 
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FIG. 1. a) Mass spectrum of dimuons produced by 800 
GeV protons incident on a copper dump. The fit described in 
the text for the Drell-Yan dimuon continuum, the T(1S), and 
the sum of the T(2S) and T(3S) resonances is also shown, b) 
The ratio of the data to the Monte Carlo-generated Drell-Yan 
continuum events; the ratio of the T(1S) and the T(2S+3S) 
generated events to the continuum fit is also shown. 

the nucleon-nucleon center-of-mass frame), and pr < 4.0 
GeV/c (pt is the transverse momentum of the dimuon). 

For this measurement the currents of the two spec- 
trometer magnets were set to 4200 A and 4265 A, their 
maximum excitation, which produced a spectrometer ac- 
ceptance that decreased rapidly for dimuon masses below 
8 GeV. Figure 1 shows the observed dimuon mass spec- 
trum from 8.1 GeV to 15.0 GeV dimuon mass. The com- 
ponents of a fit described below are also indicated. The 
smooth continuum of dimuons under the bottomonium 
peaks arises from the production of dimuons via quark- 
antiquark annihilation, the Drell-Yan process |ll]]. The 
experimentally observed width of the intrinsically nar- 
row onium states arises from muon multiple scattering 
and energy loss in the 4m-thick copper target. 

The Drell-Yan dimuon continuum is described well 
with a PQCD calculation |Q incorporating a recent 
MRST determination of the proton structure func- 
tions ]l^j . The yield of Drell-Yan dimuons is modeled 
with a Monte Carlo simulation of the apparatus that 
generates events as a function of dimuon px and the 
apparent fractional momenta, x\ and x%, of the anni- 
hilating quark-antiquark pair (where sx\x% = rn 2 and 



FIG. 2. The observed xf and pr distributions of the 
T(1S). The data distributions are formed by subtracting the 
Monte Carlo-generated Drell-Yan dimuons, and the generated 
T(2S+3S) dimuon decays, from the observed dimuon data. 
The Monte Carlo-generated T(1S) decay spectra are shown 
for comparison. A similar comparison is included for the sum 
of the T(2S) and T(3S) states 

x\ — X2 = xf] m is the dimuon mass, s is the center- 
of-mass energy squared). A standard parametrization of 
the Drell-Yan production cross section versus px was fit 
to the data |T3] . Drell-Yan virtual photons are produced 
transversely polarized and hence their dimuon decay is 
predicted to yield a 1 + cos 2 (0) angular distribution. 

Since the mass of a bottomonium state is fixed, the 
production of a bottomonium state is a function of pt 
and Xf only. The functional form of the production dis- 
tributions can be found from the data directly. Due to 
the 330 MeV rms mass resolution of this measurement, 
we cannot resolve the 2S and 3S states. It has previously 
been observed that the pt and xf distributions of the 
T(2S) and T(3S) states are very similar juj. Thus, in 
our fits to the data to extract the decay angular distri- 
butions, we assume that the 2S and 3S states have the 
same px and xp distributions. However, we note that 
the results in this paper are insensitive to this assump- 
tion within statistics. 

We generated twice as many accepted Monte Carlo 
events as were observed in the data. The Drell-Yan 
dimuon continuum was generated using PQCD with 
MRST parton distributions |D| , a shape versus pr that 
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FIG. 3. a) Decay angular distribution of T(1S) dimuon 
decays, formed by subtracting the fit contributions of the 
Drell-Yan, T(2S), and T(3S) decays from the data (in the 
bin 8.8 < m^- < 10.0 GeV and p T > 1.8 GeV/c). 
A fit to the form 1 + a cos 2 (9) is superimposed, b) The 
corresponding decay distribution for T(2S+3S) decays (for 
10.0 < m M+M - < 11.1 GeV and p T > 1.8 GeV/c). 

fit the data, and a transverse polarization of 100%. The 
Drell-Yan continuum events were then weighted with 
quadratic polynomial functions of x\ and X2 to match the 
data exactly. The weighting polynomials (which varied in 
value from 0.85 to 1.15) correct for small inaccuracies in 
the modelling of the apparatus and for variations of the 
p-Cu cross section from the PQCD prediction (there are 
known to be small nuclear effects in Drell-Yan dimuon 
yields JT^]). The weighting is important since acceptance 
correlations between muon momenta and dimuon decay 
angle could lead to a false polarization signal if the ob- 
served yield is not modelled correctly versus xf and px- 

The Monte Carlo simulation of the bottomonium 
states generated unpolarized T(1S), T(2S) and T(3S) 
events with px and xp shapes and a relative 1S/2S/3S 
weight that matched the data. In the final fit to the data, 
the polarization parameter a of both the Drell-Yan con- 
tinuum and the bottomonium resonances was allowed to 
vary. The polarizations of the T(2S) and T(3S) states 
were set equal in the fit after attempts to assign different 
polarizations to these two states led to large, negatively 
correlated statistical errors on a for the two states (con- 
sistent with the limited resolution mentioned above). 

The final fit independently varied the shapes of the pro- 



duction distributions and the polarizations of the Drell- 
Yan, T(1S) and T(2S+3S) generated events to match the 
data. The shapes of the production distributions agree, 
within errors, with those obtained earlier pj]. Figure la 
shows the results of the fit versus dimuon mass for all the 
data. Figure lb shows the ratio of the data and gener- 
ated resonances to the generated Drell-Yan events. The 
separation of the T(1S) from the combined T(2S) and 
T(3S) states is sufficient to yield a stable fit. 

In figure 2 we show the xf and px distributions ob- 
served for the T(1S) (8.8 < m M+M - < 10.0 GeV) data 
along with the fitted Monte Carlo distributions for the 
T(1S). The data spectra are obtained by subtracting the 
Monte Carlo fit distributions for the Drell-Yan, T(2S), 
and T(3S) dimuons from the data. The figure also in- 
cludes similar curves for the sum of the T(2S) and T(3S) 
states. The acceptance varies more slowly than the ob- 
served event yield versus either xp or px', the average xf 
of either the analysed T(1S) or T(2S+3S) data is 0.23 
and the average px is 1.3 GeV/c. 

Figure 3 shows the angular distributions, in one of four 
px bins, for the T(1S) decays and for the sum of the 
T(2S) and T(3S) decays. Each point in figure 3a shows 
the data in a mass bin 8.8 < m^^- < 10.0 GeV with 
the Monte Carlo-generated contributions from Drell-Yan 
dimuons, T(2S) decays, and T(3S) decays subtracted 
away. Similarly, figure 3b shows the data in a mass bin 



from 10.0 < 



< 11.1 GeV minus the Monte Carlo- 



generated Drell-Yan and T(1S) events. The expected 
1 + a cos 2 (8) decay angle distribution fits well in both 
cases. The x 2 /DF of the fits are 0.7 and 1.2 respectively. 

The values of a arising from the combined produc- 
tion distribution and decay angular distribution fit in 
the Drell-Yan sideband and two onium mass regions for 
4 bins in px (bin boundaries at px — 0.0, 0.8, 1.3, 1.8 
and 4.0 GeV/c) are shown in figure 4a. The results ver- 
sus xp (4 bins, boundaries at 0.0, 0.12, 0.23, 0.35 and 
0.6) are shown in figure 4b. The points are plotted at 
the cross-section-weighted average value of the abscissa. 
A systematic error in a of ± 0.06 should be added to 
the values of the onium polarizations in figure 4. This 
was estimated by varying the form of the fitting function 
for the Drell-Yan continuum dimuons and by varying the 
width of the mass bins used to fit the onium resonances 
and the Drell-Yan continuum. 

The observed polarization of the Drell-Yan continuum 
dimuons is consistent with 100% transverse polarization 
in all bins and with previous measurements pjj. The 
Drell-Yan sidebands have 0.2 < xi < 0.8 and 0.06 < 
X2 < 0.4, a region where no significant nuclear shadowing 
is observed |l7]] . A fit to the Drell-Yan sideband data (for 
all xf andpr) yields a — 1.008±0.016 with an estimated 
systematic error of ±0.020. 

The T(1S) data show almost no polarization at small 
xp and px- The data show a finite transverse polar- 
ization at either large px or at large xf (there are no 
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FIG. 4. a) a versus pr for the Drell-Yan sidebands 
(8.1 < m M+M - < 8.45 GeV and 11.1 < ra M+fl - < 15.0 
GeV), T(1S) (8.8 < m M+M - < 10.0 GeV), and T(2S+3S) 
(10.0 < < 11.1 GeV). b) a versus xf for the same 

mass regions. The errors shown are statistical, there is an ad- 
ditional systematic error not shown of 0.02 in a for Drell-Yan 
polarizations and 0.06 in a for onium polarizations. 

significant xf versus px production distribution correla- 
tions observed in the data). This observation disagrees 
with an NRQCD calculation that predicts a polariza- 
tion of 0.28 to 0.31 at our energies R|. If we fit the 
IS state for a polarization independent of xf and pr, we 
get a = 0.07 ± 0.04. 

The observation that the polarization of the cross- 
section-weighted average of the 2S+3S states is much 
larger than that of the IS state at all xf and px contrasts 
sharply with what is seen in the charmonium system || . 
Although an NRQCD calculation j?J predicts that feed- 
down decays from higher S, P, and D upsilon states dilute 
the polarization of the IS state, we can find no explicit 
calculation of the polarization expected for the 2S or 3S 
state. 

In the kinematic range 0.0 < xp < 0.6 and 
Pt < 4.0 GeV/c, the fit to the data yields a ratio of 
T(2S+3S)/T(1S) events of 0.50 ± 0.01. A separate 3- 
peak fit yielded an overall ratio of T(3S) to T(2S) events 
of 0.46 ± 0.03 consistent with previous high resolution 
measurements [Q. Note that even if the T(3S) were 
100% polarized, the T(2S) must be at least 35% polar- 
ized to yield the observed polarizations of the combined 



peaks. Likewise, if the T(2S) were 100% polarized, the 
T(3S) must have significant positive polarization in most 
bins. 
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